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ABSTRACT 
The  problem  of  transporting  recruits  from  Navy  Recruiting  Stations 
to  Recruit  Training  Centers  in  such  a  way  as  to  minimize  total  cost  is 
subjected  to  solution  by  standard  linear  programming  methods.  A  sensi- 
tivity analysis  of  the  variation  of  total  cost  for  different  modes  of 
transportation,  as  influenced  by  the  variation  of  the  proportion  of 
recruits  going  to  each  of  the  two  Recruit  Training  Centers,  is  carried 
out.  A  probabilistic  model  of  the  occurrence  of  a  given  number  of 
enlistments,  up  to  a  fixed  quota,  at  each  Recruiting  Station  is  formu- 
la ted  and  a  simulation  of  the  transportation  problem  arising  is  carried 
out.  A  rationale  is  developed  for  the  derivation  of  a  realistic 
estimator  for  planning,  and  some  of  the  difficulties  of  determining 
a  best  estimator  are  discussed.  A  computer  program  for  approximating 
the  expected  loss  for  any  estimator  is  presented.  Some  other  types  of 
approach  to  the  problem  are  indicated. 
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1.        Introduction. 

In  today's  Navy  considerable  emphasis  is  being  put   on  cost  re- 
ductions  of  all  forms.      One  important  example   of  an  area  in  which 
significant  cost  reductions  might  be  made  is  that  of  the   transportation 
of  recruits  from  Navy  Recruiting  Stations   (RS's)   to  Recruit  Training 
Centers    (RTC's).     The   large   number   of  such  "shipments"  made  indicate 
that  close  study  of  the  costs  are   needed „     Accordingly,   the  Bureau  of 
Naval  Personnel,   in  their  letter  PERS  B-61   KAD   of  29  May  1963,   re- 
quested investigation  of  methods  to  reduce  these   costs.     This  paper  is 
presented  as  an  illustration  of  a  practical  application  of  the  basic 
linear  programming  method   of  solving  transportation  problems  as  applied 
to  naval  recruits  from  the  time   of  their  enlistment  until  their  arrival 
at  the  appropriate  RTC. 

There  are  thirty-nine  main  Naval  Recruiting  Stations  and  nineteen 
class  "A"  substations,  a  total  of  fifty- eight  recruit  sources.     Each  RS 
is  assigned,  as  a  quota,  a  certain  percentage   of  total  number  of  recruits 
to  be  included  by  the  Navy  Recruiting  System  for  the  year.     To  receive 
the  recruits  there  are  two  RTC's  each   of  which  is  assigned  to  receive 
half  the   total  number  of  recruits  enlisted.     Recruits   can,   in  principle, 
be  sent  from  any  RS  to  either  RTC  by  any  of  three  modes   of  transporta- 
tions    air,   rail,    or  bus,  all  of  which  cost  different  amounts.     In 
addition,   there  are  differences  in  cost  which  are  a  result  of  obtaining 
group  rates.      In  order  to  analyze  the   cost  of  transportation  more   than 
superficially,    it  is   necessary  to  include  costs   other  than  those   of 
tickets,   e.g.,   meals,   miscellaneous    (limousine  service,   etc.)  and  pay 
of  recruits  during  transit.     The  desired  solution  of  the   problem 
should  give   the   optimum  number  of  recruits  to  be  sent  to  each  RTC  from 
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each  RS  and  should  determine  the  modes  by  which  they  should  -be  sente 
In  addition,  since  fare  rates,  cost  of  living,  pay  rates,  RS  quotas, 
and  RTC  quotas  may  all  vary,  it  is  necessary  to  determine  the  effects 
of  all  of  these  changes  on  the  optimum  transportation  allocations. 
Determinations  of  these  effects  are  known  as  sensitivity  analyses. 

This  problem,  with  slight  modifications,  may  be  formulated  as  the 
linear  programming  transportation  problem.  This  formulation  of  the 
problem  is  that  of  the  shipment  of  a  homogeneous  product  from  a  stated 
number  of  sources  to  a  stated  (possibly  different)  number  of  receiving 
points  (destinations) ;  the  amount  to  be  shipped  from  each  source  and 
the  amount  to  be  shipped  to  each  destination  is  predetermined;  however, 
the  amount  going  from  a  specific  source  to  a  specific  destination  is 
not.  The  cost  of  shipment  of  one  item  from  any  source  to  any  destina- 
tion is  a  necessary  input.  The  solution  determines  the  amounts  to  be 
shipped  from  each  source  to  each  destination  in  such  a  manner  as  to 
minimize  total  cost,  and  gives  the  minimum  total  cost. 

In  this  paper  we  have  made  use  of  a  previously  written  FORTRAN 
program  of  the  simplex  algorithm  to  solve  (on  a  CDC  1604.  computer)  the 
basic  transportation  problem  as  posed  by  the  Bureau  of  Naval  Personnel. 
In  addition  we  have  investigated  the  changes  in  the  solution  introduced 
by  varying  certain  of  the  constraints  (e.g.,  the  restriction  that  cj0% 
of  the  men  go  to  each  Recruit  Training  Center).  Further,  an  attempt  has 
been  made  to  come  to  grips  with  the  more  realistic  problem  in  which  the 
number  of  recruits  desired  is  specified,  but  the  number  which  will 
actually  enlist  can  be  described  only  probabilistically. 

Since  some  of  the  data  were  unavailable  at  the  time  this  paper 
was  written,  those  values  were  supplied  by  educated  guess.  Further, 


the  actual  number  of  men  in  the  quotas  were  not  known;  the  percentages 
of  the  total  were  used  instead.  Thus,  this  paper  is  a  demonstration 
of  a  method  which  becomes  a  true  solution  (within  the  limits  of  its 
assumptions)  only  when  the  true  values  are  supplied. 


2.   Problem  Formulation. 

In  this  (classical)  statement  of  the  linear  programming  transpor- 
tation problem  the  following  information  is  necessary  as  inputs 

It  is  desired  to  ship  a  homogeneous  product  from  each  of  a  given 
number  of  destinations  in  such  a  manner  that  each  destination  receives 
an  exact  desired  amount  (which  can  be  different  for  each  destination). 
Further,  this  shipment  is  to  be  carried  out  in  such  a  manner  that  each 
origin  sends  an  exact  desired  amount  which  can  be  different  for  each 
origin.  The  total  amount  that  is  desired  to  send  is  equal  to  the  total 
amount  that  is  desired  to  receive.  No  shipment  takes  place  from  destina- 
tions to  origins.  The  problem  is  to  be  solved  in  such  a  way  as  to  mini- 
mize the  total  cost  of  all  this  transportation.  In  other  problems  some 
of  the  above  equalities  may  be  inequalities,  but  in  our  problem  equal- 
ities are  adequate. 

To  solve  our  problem  realistically,  account  must  be  taken  of  the 
fact  that  seldom  will  money  savings  be  the  only  governing  criterion  in 
any  problem  involving  people.  Hence,  it  will  be  necessary  to  provide 
solutions  to  the  transportation  problem  under  consideration,  which  for 
example,  might  be  to  prevent  trip  length  from  being  excessive.  Suppose 
a  decision  were  made  (for  reasons  of  morale,  say)  that  no  recruit  should 
spend  more  than  thirty-six  hours  in  travel.  Obviously  for  some  RS's 
the  cheapest  mode  of  transportation,  which  is  also  the  slowest,  would 
no  longer  be  usable,  and  some  costs  would  be  changed.  This  would  raise 
the  total  cost  of  transporting  recruits  and  might  well  change  the  pattern 
of  shipments,  i.e.,  how  many  men  are  sent  between  which  sources  and 
destinations,  for  minimizing  total  cost  using  the  revised  costs.  Another 
variation  introduced  into  the  problem  is  that  caused  by  the  availability 


in  some  areas  of  reduced  rates  when  drafts  of  more  than  a  certain 
number  of  men  are  sent  at  one  time.  A  solution  of  the  problem  using 
the  assumption  that  advantage  is  taken  of  all  such  available  group 
rates  can  be  justified  by  assuming  that  enlistees  can  be  told  to  report 
for  transportation  a  long  enough  time  after  first  enlisting  to  ensure 
that  the  station  has  enough  recruits  to  obtain  group  rates.  It  should 
be  noted  that  various  standard  queueing  theory  models  of  this  situation 
might  be  used  to  select  appropriate  delays. 

The  input  used  in  the  problems  solved  for  this  paper  were  obtained 
as  follows; 

Costs  of  air,  train,  and  bus  tickets  from  each  RS  to  each  RTC 
(San  Diego  and  Great  Lakes)  were  requested  from  each  RS.  In  addition, 
the  RS's  were  asked  to  supply  meal  costs,  miscellaneous  costs,  and  time 
enroute  for  each  mode.  Not  all  requests  were  answered,  so  missing  values 
were  supplied  by  estimates.  In  the  data  supplied,  the  numbers  of  men 
required  at  the  destinations  and  available  at  the  origins  were  stated 
in  percents,  the  information  having  been  supplied  by  the  Bureau  of  Naval 
Personnel.  The  letter  prompting  this  study  specified  that  each  Recruit 
Training  Center  should  receive  the  same  number  of  recruits,  giving  the 
effective  determination  of  that  input  as  fifty  percent. 

Once  having  solved  the  basic  problem  for  the  various  conditions  of 
restriction  on  cost  (maximum  cost,  minimum  single  fare,  minimum  group 
fare,  minimum  with  time  limit,  etc.),  it  would  be  of  interest  to  inves- 
tigate the  effect  on  the  optimal  solution  when  changes  are  made  in  the 
system.  Since  in  the  transportation  problem,  the  addition  of  new 
variables  (RS's  or  RTC's)  requires  not  a  variation  of  the  same  problem, 
but  the  solution  of  an  entire  new  problem,  such  additions  will  not  be 


investigated.  The  changes  in  the  system  which  appear  to  lead  to  fruit- 
ful investigations  are  changes  in  the  costs,  in  the  proportion  of  re- 
cruits assigned  to  each  RTC,  and  in  the  numbers  available  at  each  RS. 
For  a  simplified  mathematical  statement  of  the  foregoing  formulation, 
see  Appendix  A. 


3.       Sensitivity  Analysis  Using  Cost  Inputs  as  Parameters . 
3.1     The  following  assumptions  are  used  for  this  analysis  % 

First  -  The  statistics  presented  in  the  tables  are  NOT  OFFICIAL 
NAVY  statistics  and  the   conclusions  from  this  analysis  are  strictly  the 
personal   opinions   of  the  authors . 

Second  -  It  was  stated  in  the  introduction   of  this  paper  that  there 
are  thirty-nine  main  RS's;   however,   this  analysis  was  conducted  using 
forty  since  we  believed  at  the  time  that  Baltimore  was  a  main  Recruiting 
Station.     This  analysis  was  completed  prior  to  the  discovery  of  this 
oversight.     This   oversight  has  no  effect  on  the   conclusions   or  the 
validity  of  the  described  method  for  solving  the   problem.     Hence,   for 
this  analysis,   recruits  are  sent  from  forty  RS  to  two  RTC's  in  accordance 
with  the  travel  rates   or  cost  coefficients  listed  in  Table   I,   page   25. 
These  rates  are  constantly  changing '9   therefore,    this  table  must  be  "up- 
dated"  to  be  applicable  to  a  particular  time.     Computed  in  these  rates 
are  subsistence  which  includes  meals,   berthing,   limousine  service,   etc., 
and  pay  during  travel  which  amounts  to  twelve  cents  per  hour. 

Third  -  The   overall  annual  Navy  recruiting  quota  is  distributed 
among  the  forty  Recruiting  Stations  as  listed  in  Table  II,   page   26. 
These  percentages  were   computed  by  multiplying  the  annual  district  quota 
by  the  individual  Recruiting  Station's  annual  quota.     For  example;, 
The  first  recruiting  district  is  responsible  for  16.52$  of  the   overall 
annual  Navy  quotau      Of  this   percent,   Boston  is  responsible  for  27.00% 
within  the  first  district;   therefore,   Boston  is  responsible  for  27.00 
times   16.52   or  4.4-6/6  of  the   overall  annual  Navy  quota.     These  percentage 
assignments  are  redistributed  slightly  each  year;   therefore,   this  table 
also  should  be  "up-dated"   to  be  applicable  to  a  particular  year. 


Furthermore,  some  of  the  RS's  have  class  "A"'  substations  which  further 
breaks  the  percentage  down;  however,,  for  the  problem  discussed  in  this 
paper  only  the  forty  main  recruiting  stations  were  considered  as  origins 
because  the  percentage  breakdowns  to  the  class  nAn   substations  was  not 
available.  Again,  this  does  not  change  the  basic  method  of  solving  the 
problem. 

Fourth  -  It  is  emphasized  most  emphatically  that  this  analysis  is 
based  only  on  costs  Incurred  during  transportation  and  ignores  any  con- 
sideration of  costs  incurred  before  arrival  of  candidates  at  RS's  (i.e., 
Recruiting  Stations  operating  costs)  or  after  arrival  at  RTC'c  .i.e., 
costs  of  operating  the  RTC's). 

The  sensitivity  analysis  was  conducted  by  using  a  fixed  set  of 
cost  input  parameters  and  varying  the  percent  recruit  input  to  RTC 
Great  Lakes  from  0%   to  100%,,  Since  the  total  of  100%  recruits  is  sent 
to  RTC  Great  Lakes  and  RTC  San  Diego,  it  is  evident  that  RTC  San  Diego 
must  receive  the  complement  of  those  recruits  sent  to  RTC  Great  Lakes. 
When  this  approach  was  used  a  particular  sequence  of  recruiting  stations 
and  corresponding  percentages  resulted 5  these  percentages  were  the  levels 
of  input  to  RTC  Great  Lakes  at  which  the  corresponding  RS's  commenced 
sending  recruits  to  RTC  Great  Lakes.   On  each  of  the  following  four  pages 
the  parameter  of  cost  input  is  defined.  Below  each  is  tabulated  the 
applicable  sequence  of  RS's. 


Max  Cost  Inputs?  Max  Cost  Inputs  are  defined  as  the  maximum 
transportation  costs,  regardless  of  mode,  between  a  RS  and  the  two 
Recruit  Training  Centers.  These  inputs  were  used  to  put  an  upper  bound 
on  the  costs.  The  use  of  these  cost  inputs  resulted  in  transporting 
the  recruits  via  all  three  modes $  however ,  each  recruit  traveled  via 
only  a  single  mode. 

Sequence  of  RS's  and  the  Corresponding  Percentage  at 
Which  Each  RS  Commences  to  Send  Recruits  to  RTC  Great  lakes 


New  York  City 

1% 

Detroit 

40% 

Raleigh 

68% 

Milwaukee 

8% 

Louisville  j, 

Ky 

44% 

Des  Moines 

70% 

Chicago 

10% 

Columbus 

46% 

New  Orleans 

72% 

Cleveland 

U% 

Jacksonville 

,Fla 

47% 

Kansas  City 

73% 

Indianapolis 

17% 

Richmond ,  Va 

50% 

Oklahoma  City 

76% 

Baltimore 

19% 

Boston 

51% 

Houston 

77% 

Albany 

20% 

Nashville 

55% 

Dallas 

79% 

Pittsburgh 

23% 

Cincinnati 

58% 

Denver 

81% 

Philadelphia 

26% 

St.  Louis 

59% 

Albuquerque 

83% 

Minneapolis 

31% 

Omaha 

62% 

Seattle 

85% 

Ashland,  Ky 

34% 

Little  Rock 

64% 

Portland 

87% 

Buffalo 

36% 

Macon,  Ga 

65% 

San  Francisco 

90% 

Washington,  DC 

38% 

Birmingham 

66% 

Los  Angeles 

94%  -  100% 

Columbia,  S.C. 

39% 

Air  Cost  Inputs;  Air  Cost  Inputs  are  defined  as  the  cost  to 
transport  all  recruits  only  by  air.  These  inputs  were  used  primarily 
to  arrive  at  a  cost  figure  when  minimum  time  enroute  is  the  primary 
consideration.  For  Chicago,  where  no  air  transportation  is  available 
to  RTC  Great  Lakes,  the  train  fare  was  used  under  the  assumption  that 
time  spend  enroute  was  less  than  by  bus. 

Sequence  of  RS's  and  the  Corresponding  Percentage  at 
Which  Each  RS  Commences  to  Send  Recruits  to  RTC  Great  Lakes 


Albany 

1% 

Indianapolis 

US 

New  Orleans 

69% 

Louisville 

3% 

Baltimore 

45% 

Omaha 

70% 

Cleveland 

4% 

Columbus 

47% 

Kansas  City 

72% 

New  York  City 

7% 

Nashville 

48% 

Little  Rock 

74% 

Boston 

15% 

Cincinnati 

50% 

Oklahoma  City 

76% 

Ashland 

19% 

Raleigh 

52% 

Dallas 

77% 

Ri  chmond 

21% 

Columbia 

53% 

Houston 

79% 

Philadelphia 

22% 

Jacksonville 

55% 

Denver 

81% 

Detroit 

27% 

St.  Louis 

57% 

Albuquerque 

83% 

Milwaukee 

31% 

Macon 

60% 

Seattle 

85% 

Washington,  DC 

33% 

Birmingham 

61% 

Portland 

87% 

Chicago 

34% 

Des  Moines 

63% 

San  Francisco 

90% 

Pittsburgh 

38% 

Minneapolis 

65% 

Los  Angeles 

94%  -  100% 

Buffalo 

42% 
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MinT  Cost  Inputs;  MinT  Cost  Inputs  are  defined  as  the  minimum 
costs  allowing  for  a  maximum  time  enroute  of  thirty -six  hours.  Thirty- 
six  hours  is  an  arbitrary  value  assumed  to  be  realistic  by  the  authors 
and  any  other  value  might  be  used  to  obtain  comparable  results.  With 
these  cost  inputs,  as  with  the  Max  cost  inputs,  the  transporting  of 
recruits  via  all  three  modes  resulted. 

Sequence  of  RS's  and  the  Corresponding  Percentage  at 
Which  Each  RS  Commences  to  Send  Recruits  to  RTC  Great  Lakes 


Ashland,  Ky 

1% 

Washington,  DC 

41% 

Little  Rock 

69% 

Boston,  Mass 

2% 

Macon 

42% 

New  Orleans 

70% 

Louisville,  Ky 

6% 

Indianapolis 

43% 

Omaha 

71% 

Ri  chmond 

s% 

Buffalo 

45% 

Kansas  City 

73% 

Philadelphia 

9% 

Chicago 

47% 

Denver 

76% 

New  York  City 

14% 

Burmi  ngham 

52% 

Oklahoma  City 

77% 

Albany 

21% 

Cincinnati 

54% 

Dallas 

79% 

Cleveland 

24% 

Jacksonville 

55% 

Houston 

81% 

Detroit 

21% 

St.  Louis 

58% 

Albuquerque 

83% 

Baltimore 

31% 

Raleigh 

61% 

Portland 

85% 

Milwaukee 

33% 

Columbia 

62% 

San  Francisco 

87% 

Columbus 

35% 

Des  Moines 

63% 

Seattle 

92% 

Nashville 

36% 

Minneapolis 

65% 

Los  Angeles 

94%  -  100% 

Pittsburgh 

38% 

• 
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Bus  Cost  Inputs;  Bus  Cost  Inputs  are  defined  as  the  costs  to 
transport  all  recruits  by  bus.  No  maximum  time  enroute  was  considered 
with  these  inputs.  We  realize  that  it  is  not  feasible  to  order  a 
recruit  to  ride  a  bus  for  seventy -eight  hours  from  Portland,  Maine  |-o 
RTC  San  Diego,  however,  these  inputs  also  correspond  to  the  minimum 
costs  available  and  thus,  put  a  lower  bound  on  the  costs. 

Sequence  of  RS's  and  the  Corresponding  Percentage  at 
Which  Each  RS  Commences  to  Send  Recruits  to  RTC  Great  Lakes 


New  York  City 

1% 

Columbus 

44% 

Birmingham 

68% 

Philadelphia 

8% 

Ashland 

46% 

Omaha 

70% 

Milwaukee 

13% 

Louisville 

47% 

New  Orleans 

72% 

Chicago 

15% 

Macon 

49% 

Kansas  City 

73% 

Pittsburgh 

19% 

Buffalo 

50% 

Oklahoma  City 

76% 

Baltimore 

22% 

Columbia 

52% 

Dallas 

77% 

Indianapolis 

24% 

St.  Louis 

54% 

Houston 

79% 

Detroit 

25% 

Raleigh 

56% 

Denver 

81% 

Cleveland 

30% 

Albany 

58% 

Albuquerque 

83% 

Minneapolis 

33% 

Des  Moines 

60% 

Portland 

85% 

Washington,  DC 

36% 

Nashville 

62% 

Seattle 

87% 

Richmond 

37% 

Jacksonville 

64% 

San  Francisco 

90% 

Cincinnati 

38% 

Little  Rock 

67% 

Los  Angeles 

94%  -  100% 

Boston 

40% 
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3.2  A  graph  of  percent  recruit  input  to  RTC  Great  Lakes  vs.   trans- 
portation  costs  per  one-hundred  recruits  is  given  in  Figure  I,   page   M+0 
This  graph  illustrates  how  the  variation  of  recruit  input  between  the 
two  RTC's  affects  transportation  costs „     The  curves   on  the  graph  were 
plotted  by  using  the   computed  data  found  in  Tables  III  and  IV,   pages 
27  and  28  respectively. 

3.3  A  FORTRAN  listing  of  the   linear  program  used  for  the  previous 
sensitivity  analysis  is  given  in  Appendix  E,     This  listing  is  included 
as  an  aid  for  those  who  desire   to  pursue  this  transportation  problem 
further;  a  detailed  explanation  of  this  program  is  not  considered 
essential  to  the  results   of  this  paper.     A  format  of  the   computer 
print-out  of  the  above  program  for  a  5%  interval,   starting  at  30% 
recruit  input  to  RTC  Great  Lakes,  and  using  1$  increments  and  Air 
cost  inputs  is  given  in  Tables  VI -A  through  VI-Fy   pages   30  to  35. 
Table  V,   page   29,   presents  a  number  code  which  correlates  the   P-vectors 
in  Tables  VI -A  through  VI -F,   pages   30  to  359   to  their  respective  RS's 
and  RTC's. 
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PERCENT  RECRUIT  INPUT  TO  GREAT  LAKES 
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U.       Statistical  Variation. 

In  the  previous  formulation  of  our  problem,  the  basic  assumption 
was  that  the  inputs  to  the  problem  were  known  perfectly,  but  it  is 
quite  likely  that  this  will  not  be  the  case.  Indeed,  it  is  apparent 
that  there  is  a  positive  probability  that  the  number  of  men  enlisted 
at  any  RS  will  not  be  equal  to  that  RS's  assigned  quota.  In  other  words, 
the  number  of  men  to  be  transported  from  any  RS  may  be  regarded  as  a 
random  variable.  Regarding  this  random  variable,/!),  as  a  parameter 
to  the  problem,  we  note  that  Garvin  3   points  out  that  for  any 
parameter,  B  ,  whose  probability  density  (or  mass)  function  is  known, 
the  minimum  of  the  objective  function,  z  ,  obeys  the  relation 


0  (B)]  -  zQ  (e[b])  U.1) 


where  s 


E 


B 


is  the  vector  of  expected  values  of  the  rVs. 
zQ  (B)  =  min  z(b),  the  minimum  objective  function  expressed  as  a 
function  of  B. 
The  lack  of  equality  above  makes  the  problem  of  saying  something 
meaningful  about  z  more  difficult. 

Even  if  we  could  easily  find  the  value  in  the  left  side  of  equa- 
tion (4-.1),  we  would  have  an  expression  only  for  an  ideal  situation 
which  cannot  exist  in  the  real  world.  In  fact  what  we  wish  to  know  is 
what  policy  we  should  choose  in  order  to  approach  this  ideal  more 
closely  than  by  the  use  of  any  other  policy.   One  approach  to  such 
problems  is  to  find  an  estimator  for  the  parameter  B,  and  then  to  use 
this  estimator,  B,  to  define  our  policy. 
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In  attempting  to  find  the    optimum  estimator,  B,  we  have  made 
certain  assumptions  concerning  the  uncertain   (i.ec,   statistically 
varying)   value   of  B,  which  is  a  vector  representing  the   number  of  en- 
listments at  each  RS  plus  the   number  of  enlistees  going  to  one   of  the 
RTC'So     These  assumptions  and  the  models  arising  therefrom  are  dis- 
cussed fully  in  Appendix  B. 

It  was  found   (see  Appendix  B)   that  there  is   no  simple  analytical 
solution  to  the  problem  of  determining  a  best  estimator.     The  problem 
of  even  evaluating  any  given  estimator  by  analytical  means  is  so  great 
as  to  be  impractical   to  solve;   so  the  authors  have   taken  the  approach 
of  using  the  same  assumptions  in  devising  a  Monte  Carlo  simulation. 
As  discussed  more  fully  in  Appendix  B,   it  is  assumed  that  the  number 
of  people  arriving  at  a  recruiting  station  follow  the   Poisson  probabil- 
ity law  and  that  the   number  of  those  who  actually  enlist  is  governed 
by  the  binomial  probability  law.      In  our  simulation,    this  random  amount 
is  determined  for  each  RS  and   the  problem  solved  to  find  the  least  cost 
for  that  set  of  inputs.     The  program  is  designed  to  solve  any  desired 
number  of  such  subproblems  in   order  to  find  an  approximation  to  the 
value   of  the  left  side   of  equation   (4..1). 

The   next  assumption  we  make   is  that  there   is,   conceptually,   some 
loss  in  not  solving  the   sub-problems  this  way,  and  that  by  following 
a  fixed  policy  we  will  have  to  incur  this  loss  at  least  part  of  the 
time.     By  the  use   of  an  estimator  we  can  find  a  policy „     The  policy 
is   the  transportation  pattern  developed  by  solution  to  a  sub-problem 
in  which  the  estimator  that  we   choose  is  used  as  the  input.     When  this 
policy  is  applied  to  the  sample   of  sub-problems  we   obtain,   for  each 
sub-problem,   a  difference   in  cost  which  will  be  either  zero  or  more 
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and  which  is  the   loss  for  that  particular  sub-problem.      In  our  Monte 
Carlo  simulation  we  have  evaluated   only  one  estimator  j,   the  average  B, 
Of  course,   given  any  other  estimator,  we  could  evaluate  it  also. 
Having  found  the   losses  we  can  find  the  average   losss  which  is  termed 
sample  risk.     We  can  also  find  the  sample  deviation  of  loss#   thereby 
gaining  some  idea  of  how  much  any  particular  loss  may  vary  from  this 
average.     With  larger  and  larger  numbers   of  sub-problems  we   can  more 
closely  approach  the  "true11  values  of  these  statistics.     The  results 
of  the  Monte  Carlo  simulation  are  discussed  in  section  5.3. 

Another  possible  estimator  which  might  be  used  under  certain  circum- 
stances is  discussed  in  Appendix  D. 
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5.       Results  and  Discussion. 

5.1      It  is  interesting  to  observe   that  the  sequence   of  the  RS's 
changed  when  the   cost  inputs  were   changed  for  some  reason,   such  as 
changing  the  mode   of  travel   (air  to  bus),    or  changing  the   criterion 
for  choosing  a  particular  mode   of  travel   (limited  time  enroute). 

Using  Max  Cost  Inputs,   it  is   noted  that  when  the  input  to  RTC 
Great  Lakes  is  changed  to  1%,   New  York  is  the  first  RS  to  send  recruits 
to  RTC  Great  Lakes,   and  it  sends  1%  of  its  quota.     New  York's  total 
quota  is  7.35%;   therefore,   as   the  input  to  RTC  Great  Lakes  is  increased 
by  1$  increments,   New  York  increases  its     quota  by  \%  increments  until 
its  total  quota  is  sent  there.      One  explanation  for  this   could  be   that 
New  York  is  responsible  for  the  largest  overall  Navy  quota.     However, 
referring  to  Table  I,   page  25,   it  is   observed  that  the  DIFFERENCE  in 
Max  costs  for  New  York  is  $211.02  -  $44.96  =  $166,06.     This  difference 
is   larger  than  for  any   other  city^   therefore,   it  would  seem  that  this 
combination  would  give  the  greater  savings. 

As  the  recruit  input  to  RTC  Great  Lakes  continues  to  change  by  1% 
increments  it  is   noted  that  when  the  input  reaches  8%,  Milwaukee  is  the 
next  RS  to  change  its  distribution  and  commence  sending  recruits  to 
RTC  Great  Lakes.     Why  should  Milwaukee   send  recruits  to  RTC  Great  Lakes 
before  Chicago?     Once  again  it  is  observed  that  the  DIFFERENCE  in  Max 
costs  for  Milwaukee   is  $152.99  -  $10.09  =  $U2.90,   compared  to  $U3.39 
-  $1.37  =  $142.02,   for  Chicago.     It  is  interesting  to  note  that  the 
input  to  RTC  Great  Lakes  must  reach  10%  before  Chicago  starts  sending 
recruits  to  Great  Lakes.     Although  Chicago  is  the   closest  recruiting 
station  to  Great  Lakes,    this  does  not  imply  that  the  minimum  distance 
traveled  is  going  to  give  maximum  savings. 
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Again  referring  to  Table  I,  page  25,  only  this  time  using  the 
Air  or  Bus  costs,  it  is  observed  that  the  first  RS  to  commence  sending 
recruits  to  both  RTC's  is  not  New  York,  but  that  city  that  has  the 
greatest  DIFFERENCE  in  the  appropriate  costs.  When  that  RS  has  com- 
menced sending  its  total  quota  to  RTC  Great  Lakes,  the  next  city  with 
the  greatest  DIFFERENCE  in  costs  commences  to  send  recruits  to  both 
RTC's.  This  method  continues  until  all  RS's  are  sending  their  total 
quota  to  RTC  Great  Lakes.  Therefore,  when  the  problem  is  linear  and 
involves  two  destinations,  as  in  this  problem,  we  see  that  it  is  the 
DIFFERENCE  in  costs  of  the  applicable  mode  of  travel  to  the  two  des- 
tinations that  is  the  controlling  factor. 

5.2  The  graph  in  Figure  1,  page  14-,  illustrates  that  there  is  a  de- 
crease in  travel  costs  up  to  an  82$  input  of  recruits  to  RTC  Great 
Lakes  and  then  the  costs  begin  to  increase.  In  the  case  of  Max  and 
MinT  cost  inputs,  where  it  is  recalled  that  mixed  travel  modes  were 
used,  the  minimum  occurred  at  roughly  the  same  percentage.  Thus,  the 
decrease  occurs  regardless  of  the  mode  of  transportation. 

Why  should  the  minimum  costs  occur  at  such  a  high  percentage 
input  to  RTC  Great  Lakes?  It  can  be  conjectured  that  one  possible 
explanation  might  be? 

1 .  It  costs  less  to  transport  recruits  from  all  parts  of  the 
country  to  RTC  Great  Lakes  than  to  RTC  San  Diego  because  RTC 
Great  Lakes  is  geographically  more  centrally  located. 

2.  The  more  populous  areas  are  responsible  for  a  larger  overall 
annual  recruiting  quota  and  the  majority  of  these  areas  are 
located  nearer  RTC  Great  Lakes  than  RTC  San  Diego. 
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We  would  like  to  emphasize  at  this  point  that  we  have  discussed 
only  that  the   transportation  costs  to  RTC  Great  Lakes  are   less  than 
those  to  RTC  San  Diego  up  to  a  certain  percentage   of  recruit  input. 
We  have  made   no  cost  analyses   of  sending  recruits  through  the  two 
Recruit  Training  Centers. 

5.3     In  testing  the  Monte  Carlo  approach  to  this  problem,   first  a 
sample   of  fifty  randomly  generated  B-vectors  were  analyzed.     The  inputs 
to  the  random  vector  generator  sub-program  were, 

1)  RS  quotas  expressed  as  integer  numbers, 

2)  the  average   number  of  interviews  held  per  day  at  each  RS,   and 

3)  the   individual  probability  of  success    (obtaining  an  enlistment) 
per  interview  for  each  RS. 

Values  of  the   last  two  parameters  were  quite  arbitrarily  chosen,   the 
average  interview  number  being  made  approximately  a  fourth  more  than 
the  annual  quota  divided  by  365.     The  probability  of  success  at  an 
interview  was   chosen  arbitrarily  over  a  range  from  55%  to  80%.     A  table 
of  these  inputs  is  given  in  Table  VII-B.     The  vectors  resulting  had 
entries   of  a  fair  number  of  stations   (about  half  in  most  cases)  which 
fulfilled  their  quota.     A  typical  vector  resulting  is  shown  in  Table  VII-C. 
The   other  per  cent-of -quotas-enlisted  ranged  down  to  60>t  but  the   total 
quota  was  filled  to  a  level  never  less   than  93%  and  never  as  great  as 
95%  in  the  fifty  samples.     Further,   the  solution  pattern  of  transporta- 
tion was  the  same  for  each  random  vector,   for  these  arbitrary  parameter 
values.     Typical  results  of  the  solution  of  one  vector  are  given  in 
Table  VII-D.     In  short,   the  solution  was  quite   insensitive   to  changes 
generated  in  a  small  sample  by  these  parameter  values,  and  the  sample 
deviation  of  the  random  variable  p   for  the  RTC  was  small.     Because 
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of  this  result,  it  was  decided  that  more  runs  would  be  desirable  to 
determine  whether  different  parameter  values  might  lead  to  very  different 
results.  Accordingly,  four  samples  were  computed,  with  fifteen  random 
vectors  each  and  with  combinations  of  parameters  that  differed  greatly. 
The  average  number  of  interview  for  each  RS  was  given  the  value  of  the 
quota  of  that  RS  divided  by  365  for  two  runs  and  twice  that  much  for 
two  others.  The  probability  of  interview  success  was  fixed  by  55%  for 
all  stations  for  two  runs  and  65%  for  the  other  two  runs.  These  were, 
of  course,  combined  in  such  a  way  as  to  get  all  four  combinations. 

Finally,  our  simulation  also  produced  sample  mean  loss,  an  approx- 
imation to  the  risk  for  the  given  estimator,  and  the  sample  deviation 
of  that  loss.  Examples  of  these  results  are  included  in  Table  VIX-D. 
When  low  averages  and  low  probabilities  (55%)  were  used,  the  percent 
of  total  quota  filled  never  dropped  below  56%  and  never  reached  57%. 
Individual  percent  of  quotas  ranged  from  39%  to  73%.  When  low  averages 
and  high  probabilities  (75%)  were  used  the  percent  of  total  quota 
filled  ranged  between  75%  and  77%  while  individual  quotas  were  filled 
to  levels  from  55%  to  96%.  For  the  combination  of  high  average  number 
of  arrivals  and  low  probability  of  enlistments  total  quota  fulfillment 
never  fell  below  99%  while  individual  station  quotas  were  filled  to 
levels  as  low  as  80%  for  the  very  few  in  each  sample  not  fulfilling 
their  quotas.  The  combination  of  high  probability  and  high  average 
gave  100%  fulfillment  for  all  fifteen  samples. 

The  lack  of  variation  in  the  percent  of  total  quota  filled  can 
be  viewed  as  a  substantiation  of  the  validity  of  the  statement  that 
for  large  samples  /jm+^  can  be  approximated  by  a  normal  random  variable 
independent  of  the  other  m  random  variables. 
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This  lack  of  variation  had  one  other  notable  effect.  The  mean 
B  vector,  B,  proved  to  be  a  very  good  estimator j  there  were,  in  fact, no 
samples  on  which  the  result  was  non-optimal.  It  should  be  noted  that 
should  the  use  of  the  transformation  P   result  in  a  sample  with  a  non- 
optimal  solution,  our  program  will  compute  the  loss  on  that  sample. 
Therefore,  not  only  was  B  a  good  estimator  for  our  selected  inputs,  but 
it  is  also  obvious  that  for  these  inputs ,  a  static  model  is  adequate. 
Of  course,  use  of  different  quotas  and  different  parameters  might  cause 
the  probability  mass  functions  governing  the  inputs  to  "spread"  so  much 
that  non-optimal  solutions  might  result  in  enough  cases  to  make  a  dyna- 
mic model  desirable. 

5.4-  In  the  probability  model  developed  herein  many  simplifying  assump- 
tions were  made  which  might  have  an  effect  on  the  conclusions.  Some 
of  them  are  discussed  below. 

1.  It  was  assumed  that  quotas  assigned  to  RS's  were  fixed.  It 
has  been  learned  that  such  quotas  may  in  fact  be  changed 
during  the  year  in  order  to  capitalize  on  peculiar  seasonal 
or  regional  situations  effecting  recruitment. 

2.  It  was  assumed  that  the  number  of  interviews  held  by  a  recruit- 
ing station  could  be  described  by  a  Poisson  probability  mass 
function.  In  using  the  Poisson  distribution,  one  assumes  that 
the  probability  of  more  than  one  event  occurring  in  time  At 
is  of  smaller  order  than  At  so  that  as  At  approaches  zero 
the  probability  of  two  or  more  events  occurring  in  At  also 

is  negligible.  In  practice  recruiters  often  give  talks  to 
large  groups  of  high  school  or  college  students.  This  might 
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be  said  to  be  more  than  one  interview;  however,  it  is  believed 
that  the  prospective  enlistees  must  still  enter  a  queue  of 
some  sort  before  making  a  commitment  to  enlist,  so  that  the 
mass  lecture  may  be  regarded  as  one  of  the  factors  determining 
0( ,  the  average  number  of  interviews  for  the  j™1  RS  rather 
than  as  an  actual  multiple  interview. 

3.  Probabilistic  independence  of  the  occurrence  of  interviews 
was  assumed.  Arguments  contrary  to  this  assumption  have  been 
acknowledged  above  and  some  defense  given.  As  a  further  de- 
fense, we  note  that  the  amount  of  personal  effect  required  of 
the  recruiter  is  a  factor  tending  to  promote  the  statistical 
independence  of  these  random  variables. 

4.  It  has  been  assumed  that  the  parameters  of  the  Poisson  and 
binomial  distributions  used  could  be  determined  for  each  RS. 
They  cannot  be  completely  determined,  of  course,  but  it  is 
likely  that  enough  data  can  be  obtained  from  RS  records  to 
arrive  at  very  good  estimates. 

Although  some  of  these  objections  are  forceful,  none  is  believed  to  be 
strong  enough  to  invalidate  the  usefulness  of  the  results. 
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6.   Conclusions. 

It  can  be  concluded  from  this  investigation  that: 

1 .  We  have  been  able  to  devise  a  method  enabling  us  to  simulate 
(at  least  approximately)  the  uncertainties  caused  by  the 
probabilistic  nature  of  numbers  of  enlistments  and  thus,  can 
consider  our  model  to  be  more  than  a  mere  laboratory  affair. 

2.  We  have  been  able  to  develop  a  rationale  for  the  derivation 
of  an  estimator  and  have  been  able  to  obtain  an  approximation 
of  the  success  we  may  expect  in  solving  our  problem  relative 
to  an  ideal  solution,  given  any  particular  estimator. 

3.  By  means  of  simulation,  we  have  found  that  the  mean  input 
vector,  B,  is  an  estimator  which  yields  optimal  solutions  for 
a  very  large  number  of  samples  and  for  widely  varying  para- 
meters. We  have  seen  that  this  indicates  the  adequacy  of  a 
static  model. 

A.  We  have  indicated  another  promising  estimator  that  might  be 
developed  that  could  have  considerable  appeal  to  a  planner 
under  certain  circumstances . 

5.  We  have  been  able  to  devise  a  method  allowing  us  to  present 
to  the  planner  an  assessment  of  the  cost  of  certain  decisions, 
namely  the  choice  of  assigning  various  quotas  to  each  of  the 
two  RTC's,  the  choice  of  putting  a  limit  on  travel  time,  and 
the  choice  of  an  estimator. 

6.  We  have  found  the  simple  principle  underlying  the  choice  of 
routes  to  be  used  by  the  RS's,  that  of  maximum  cost  DIFFERENCE, 
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TABLE     III 

TABULATION    OF -MAXIMUM    TRAVEL    COSTS    PER    RECRUIT    AND 
CORRESPONDING    PERCENTAGE     INPUT. TO    GREAT    LAKES 


00 

128.  6  9 

16 

91  .04 

52 

63.  18 

78 

43.76 

01 

12  7.  0  3 

27 

39.76 

55 

62.28 

79 

43.45 

02 

125.  57 

28.  - 

88.49 

.  54 

61  .38 

30 

45.31 

0.3 

-■  2  3.  7  1 

29 

87.22 

5o 

60.4  8 

81 

43.16 

04 

122. 0  5 

30 

85.95 

56 

59.58 

82 

43.0  5 

05 

120.39 

31 

84.75 

57 

53.69 

83 

43.01 

06 

118.73 

:  32 

3  5.59 

58 

57.79 

84 

45.22 

07 

117.0  7 

33 

82.  4  3 

59 

56.92 

85 

43.54 

OB 

115.56 

34 

3  1.2.8 

60 

56. 11 

86 

43.88 

09 

114.  13 

35 

80.  16 

61 

55.  31 

87 

44.23 

10 

112.  70 

36 

79.07 

62 

54.50 

'  88 

44.69 

1  1 

111.2  8 

37 

77.99 

65 

53.70 

89 

45.  15 

12 

109.  36 

38 

7o.90 

64 

52.99 

90 

46.  15 

13 

108.  44 

39 

75.85 

65 

52.11 

91 

47.16 

14 

'0  7.  0  2 

40 

7U.81 

66 

51  .33 

92 

48.  16 

15 

105.61 

Ul 

73.81 

67 

50.57 

9  3 

49.  17 

16 

104.  19 

42 

72.81 

68 

49.81 

94 

50.23 

17 

102.78 

43 

71  .80 

69 

49.06 

95 

51  ..51 

18 

101.4  1 

44 

70.8  0 

70 

•  48.31 

96 

52.80 

19 

100.  06 

45 

69.81 

71 

47.57 

97 

54.09 

20 

98.  7  2 

46 

68.83 

72 

46.84 

98 

55.37 

21 

97.  43 

47 

67.87 

73 

46.13 

99 

56.66 

22 

96.  14 

48 

66.90 

74 

45.45 

100 

57.95 

23 

94.  86 

49 

65.93 

75 

44.74 

24 

93.  58 

50 

65.00 

76 

44.39 

25 

92.  31 

51 

64.03 

77 

44  07 

TABULATION  OF 

AIR  TRAVEL 

COSTS 

PER  RECRUIT  AND 

SORRE 

SPONDING 

PERCENTAGE' 

INPUT 

TO  GREAT 

LAKES 

00 

102.33 

26 

76.00 

52 

51.31 

78 

34.  IS 

01 

101.  74 

27 

7  5.  00 

53 

50. 4o 

79 

34.04 

02 

100.  64 

2  3 

74.00 

54 

49.62 

80 

33.94 

03 

99.  54 

29 

72.99 

55 

48.79 

81 

33.84 

04 

98.  45 

20 

71  .99 

56 

47.93 

82 

33.82 

05 

97.  42 

31 

70.99 

57 

47.18 

8.3 

33.90 

06 

9  6.38 

32 

70.02 

58 

46.39 

84 

34.20 

07 

9  5.  34 

53 

69.05 

59 

45.60 

85 

34.62 

08 

9  4.  32 

34 

68.09 

60 

44.31 

36 

35.05 

09 

9  3.  2  9 

3  5 

67.  14 

61 

44  .03 

87 

35.49 

10 

92.  2  7 

36 

66.  13 

62 

45.27 

83 

35.95 

11 

91.  2i* 

37 

65.22 

63 

42.51   • 

39 

36.40 

12 

90.  2  1 

33 

64.26 

64 

41  .76 

90 

37.  1  1 

13 

89.  19 

39 

62.31 

65 

4  1  .02 

91 

37.81 

14 

88.  16 

40  . 

62.36 

66 

40.31 

92 

38.51 

15 

87.  1  n 

41 

61  .41 

67 

39.61 

93 

39.22 

16 

8  6.  12 

.  42 

60.46 

68 

38.90 

94 

39.96 

17 

8  5.  10 

4  3 

59.52 

69 

38.29 

95 

40.83 

13 

84.  08 

44 

58.57 

70 

57.71 

96 

41  .71 

19 

3  3.  0  6 

45 

57.63 

71 

37.13 

97 

42.59 

20 

82.  05 

^6  •" 

»?  56.6  9 

72 

36.57 

98 

43.46 

21 

81.  04 

'47 

55.76 

73 

36.02 

99 

44.34 

22 

80.  03 

48 

54.34 

74 

35.49 

100 

45.22 

23 

79.  02 

49 

55.92 

75 

55.01 

2U 

78.02 

50 

53.02 

76 

34.67 

25 

77.01 

51 

52.16 

77 

34.40 
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TABLE    IV 

rABULATION    OF    MINT    TRAVEL    COSTS    PER    RECRUIT    AND 
CORRESPONDING    PERCENTAGE     INPUT    TO    GREAT    LAKES 


00 

99.  4  3 

26 

70.05 

52 

4?  .69 

73 

24.10 

01 

98.22 

27 

68.96 

53 

42.73 

79 

23.92 

0  2 

97.  0  1 

28  ■ 

67.  fa 

'   54 

41  .73 

80 

23.73 

03 

95.  81 

29 

66.80 

55 

40.34 

8  1 

2  3.67 

04 

9U.  61 

30 

65.72 

56 

39.92 

82 

23.57 

C5 

93.  4  0 

51 

64.64 

57 

39.00 

85 

23.54 

0  6 

92.  20 

32 

65.59 

58 

38 .10  i 

84 

23.65 

07 

9  1.02 

33' 

62.54 

59 

37.20 

85 

24.03 

0* 

8  9.  8  9 

34 

61  .49 

60 

35.31  J 

86 

24.55 

09 

83.  77 

35 

60. 4u 

61 

35.46 

37 

25.24 

10 

87.66 

36 

59.4  1 

62 

34.61 

88 

26.00 

11 

8  6.  5  5 

37 

58.38 

63 

33.76 

89 

26.77 

•12 

85.  Uii 

38 

57.36 

64 

52.9.1 

90 

27.55 

13 

84.  34 

39 

56.34 

65 

32.09 

91 

28.30 

14 

8  3.  2  3 

40 

55.32 

66 

31  .27 

92 

29.  13 

15 

82.  13 

41 

54.31 

67 

60.45 

93 

29.95 

16 

81.03 

42 

53.33 

68 

29.6  3 

94 

30.31 

17 

79.  92 

4  3 

52.35 

69 

28.90 

95 

31  .75 

18 

78.  8  2 

44 

51  .38 

70 

23.2  0 

96 

32.70 

19 

77.72 

45 

50.41 

•  71 

27.52 

97 

33.64 

20 

76.6  1 

46 

4  9.45 

72 

26  . 8  5 

98 

34.58 

21 

75.  5  1 

47 

48 .  4  8 

'  7  3 

26.21 

99 

35.53 

22 

74.  42 

48 

47.52 

74 

25.60 

100 

36.47 

23 

73.  3  2 

49 

46.5  7 

75 

2M.99 

24 

72.23 

50 

45.61 

76 

24.67 

25 

71.  14 

51 

44.65 

77 

24.37 

TABULATION  OF 

BUS  TRAVEL 

COSTS 

PER  RECRUIT  AND 

:CRRESPONDING 

PERCENTAGE 

INPUT 

TO  GREAT 

LAKES 

00 

68.7  5 

26 

4^.60 

52 

33.29 

73 

22.87 

CI 

67.  8  8 

27 

47.92 

55 

32.79 

79 

22.74 

02 

67.  00 

28 

47.24 

54 

32.31 

80 

22.62 

03 

66.  12 

29 

46.56 

55 

31.33 

81 

22.51 

04 

6  5.25 

30 

45.90 

56 

31  .36 

82 

22.1*8 

05 

64.  57 

31 

45.2  5 

57 

30.39 

8  3 

22.50 

06 

63.  h9 

32 

44.59 

58 

30.43 

34 

22.61 

07 

62.  6  2 

33 

4  3.97 

59 

29.97- 

85 

22.9  5 

08 

61.  79 

34 

4i.  35 

60 

29.52 

86 

23.32 

09 

61.  0  0 

35 

42.73 

61 

2  9.06 

87 

23.69 

10 

60.  2  1 

36 

42.  1  1 

62 

2:~  .6  1 

83 

24.07 

1  1 

59.4  1 

37 

41.50 

65 

2r<.17 

89 

24.44 

12 

58.6  2 

38 

40.91 

64 

27.75 

90 

24.92 

13 

57.  8  6 

39 

40.34 

65 

27.29 

,   91 

25.41 

14 

57.  1  1 

40 

39.73 

66 

26.86 

•   92 

25.90 

15 

56.  33 

41 

39.22 

67 

'•26.4  6 

9  3 

26.38 

16 

5  5.  6  5 

42 

38.66 

68 

26.07 

94 

26.91 

17 

54.93 

43 

38.  10 

69 

25.70 

95 

27.59 

13 

54.  20 

44 

37.5  4 

70 

25.32 

96 

23.27 

19 

53.48 

45 

36.99 

71 

24.94 

97 

28.95 

20 

5  2.  7  6 

46 

36.44 

72 

24.57 

98 

29.62 

21 

52.  05 

hi 

35.90 

73 

24.21 

99 

30.30 

22 

51.  34 

48 

35.  36 

74 

23.85 

100 

30.98 

23 

50.  64 

49 

34.34 

75 

23.49 

24 

49.  96 

50 

34.32 

76 

23.21 

25 

49.  2e 

51 

3  3.80 

77 

23.01 
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TABLE  V 


NUMBER  CODE  OF  RECRUITING  STATIONS  AND  RECRUIT  TRAINING  CENTERS 


RECRJITING  STATION  GLAKES  SDIEGO 

ALBANY  PI  P2 

BOSTON  Pi  Pk 

BUFFALO  P5  P6 

•  NEW  YORK  P7  P8 

ASHLAND  P9  PlO 

BALTIMORE  Pll  P12 

LOUISVILLE  P13  P1U 

RICHMOND  PIS  P16 

WASHINGTON  P17  P18 

BIRMINGHAM  PI  9  P20 

COLUMBIA  P27      .    P22 

JACKSONVILLE  P23  P2U 

MACON  P25  P26 

NASHVILLE  P27  P2S 

RALEIGH  P29  P5C 

CINCINNATI  P31  P52 

CLEVELAND  P33  P5U 

COLUMBUS  P35  P36 

PHILADELPHIA  P37  P3S 

PITTSBURGH  P39  P40 

CHICAGO  PUi  P42 

DETROIT  PU3  PUU 

INDIANAPOLIS  PU5  pi+6 

MILWAUKEE  PU7  P48 

ST  LOUIS  PU9  P50 

DENVER  P51  P52 

DES  MOINES  P53  P5U 

KANSAS  CITY  P55  P56 

MINNEAPOLIS  P57  P58 

OMAHA  P59  P60 

ALBUQUERQUE  P6l  P62 

DALLAS  P63     ,  P6U 

HOUSTON  P65    '  P66 

LITTLE  ROCK  P67  P68 

NEW  ORLEANS  P69  P70 

OKLAHOMA  CITY  P71  -  P72 

LOS  ANGELES  '  P73  P7U 

PORTLAND  P75  P76 

SAN  FRANCISCO  P77  P78 

SEATTLE  P79  P80 
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TABLE  VI-A 


PRDBLEM  RECRUIT   «SAMPLE« 


MINIMUM  COST  OF  OBJECTIVE  FUNCTION  IS    .719862E+04 


GREAT  LAKES  PERCENT 


.300E+02 


BASIS  VECTORS  AND  COEFFICIENTS 


VECTDR 


COEFFICIENT     ( X-ZERO    COMPONENT) 


P( 
P( 
P( 
P( 
P( 


1) 
3) 
6) 
7) 
9) 


P<43) 

P(  13) 

P(  15) 

P(  IB) 

P<37) 

(22) 

(24) 

(2  6) 

(28) 

(30) 

P(32) 

P(20) 

(36) 

(12) 

(40) 

(42) 

(44) 

(46) 

P(48) 

P(50) 

(52) 

(54) 

(5  6) 

(58) 

(60) 

(62) 

(64) 

(66) 

(68) 

(70) 

(72) 

(74) 

(76) 

(78) 

(80) 

(33) 


.256060 
.446040 
.21U760 
.735140 
.  136107 
.365571 
. 139941 
.  139941 
.830700 
.488250 
.147832 
.244575 
. 135875 
.205443 
.147832 
.167400 
.205443 
. 139500 
.139941 
.306900 
.456274 
.682690 
. 166334 
.192276 
.277732 
. 180803 
. 186418 
.250429 
.338023 
. 167327 
. 143260 
.188708 
.253916 
. 136344 
.124488 
.  141284 
.620100 
.214650 
.480180 
.275070 
.292950 


E  +  0 
E  +  0 
E  +  0 
E+0 
E  +  0 
E  +  0 
E  +  0 
E+0 
E  +  00 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E+0 
E  +  0 
E  +  00 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E+0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
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TABLE  VI -B 


PR3BLEM    RECRUIT       *SAMPLE« 


MINIMUM    COST    OF    OBJECTIVE    FUNCTION    IS         .709928E+04 


GREAT    LAKES    PERCENT 


.310E+02 


BASIS  VECTORS  AND  COEFFICIENTS 


VECTDR 


COEFFICIENT  (X-ZERO  COMPONENT) 


P( 
P( 
P( 
P( 
P( 
P( 


1) 
3) 
6) 
7) 
9) 
12) 


(13) 
(15) 
(47) 
(33) 
(22) 
(37) 
(26) 
(20) 
(30) 
(32) 
(24) 
(36) 
(28) 
(40) 
(42) 
P(  18) 
P(46) 
148) 
(50) 
(52) 
(54) 
(56) 
(58) 
(60) 
(62) 
(64) 
(66) 
(68) 
(70) 
(72) 
(74) 
(76) 
(78) 
(80) 
(43) 


.2560 
.4460 
.2147 
.7351 
.1361 
.  1399 
.  1399 
.  1399 
.3173 
.2929 
.  1478 
.4882 
.1358 
.2054 
.  1478 
.1674 
.2445 
.  1395 
.205U 
.3069 
.4562 
.8307 
.  1663 
.  1605 
.2777 
.  1808 
.1864 
.2504 
.3380 
.  1673 
.  1432 
.  1887 
.2559 
.1363 
•  1244 
.1412 
.6201 
.2146 
.4801 
.2750 
.4338 


60E+0 
40E  +  0 
60E+0 
40E  +  0 
07E  +  0 
41E  +  0 
41E  +  0 
41E  +  0 
10E+00 
50E+0 
32E  +  0 
50E  +  0 
75E  +  0 
43E+0 
32E  +  0 
OOE  +  0 
75E+0 
OOE  +  0 
43E  +  0 
OOE  +  0 
74E+0 
OOE+00 
34E  +  0 
45E  +  0 
32E  +  0 
03E  +  0 
18E  +  0 
29E  +  0 
23E+0 
27E+0 
60E  +  0 
08E  +  0 
16E  +  0 
44E  +  0 
88E  +  0 
84E  +  0 
OOE+0 
50E+0 
80E  +  0 
70E  +  0 
40E+0 
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TABLE  VI -C 


PROBLEM.  RECRUIT       *SAMPLE» 


MINIMUM   COST    OF    OBJECTIVE    FUNCTION    IS         .700227E  +  C4; 


GREAT    LAKES    PERCENT 


.320E+02 


BASIS  VECTORS  AND  COEFFICIENTS 


VECTOR 


COEFFICIENT  (X-ZERO  COMPONENT) 


(  1) 
(  3) 
(47) 

(  7) 
(  9) 

(12) 
(  13) 
(15) 
(18) 
(33) 
P(22) 
P(24) 
(37) 
(20) 
(30) 
(32) 
(26) 
(36) 
P(28) 
P(UO) 
P(42) 
(  6) 
(46) 
(48) 
(50) 
(52) 
(5U) 
(56) 
(58) 
(60) 
(62) 
(6U) 
(66) 
(68) 
(70) 
(72) 
(74) 
(76) 
P(78) 
P(80) 
P(43) 


.256060E+0 

.446040E+0 

.  131731E  +  0 

.735140E+0 

. 136107E+0 

. 139941E+0 

-.139941E+0 

. 139941E+0 

.830700E+00 

.292950E+0 

.147832E+0 

.244575E+0 

.488250E+0 

.205443E+0 

.147832E+0 

.  167400E  +  0 

. 155e75E+0 

. 159500E+0 

.205443E+0 

.306900E+0 

.456274E+0 

.214760E+0 

. 166334E+0 

.605450E+00 

.277732E+0 

. 180803E+0 

.186418E+0 

.250429E+0 

.338023E+0 

. 167327E+0 

.143260E+0 

. 138708E+0 

.253916E+0 

. 136344E+0 

.124488E+0 

.  141284E  +  0 

.620100E+0 

.21465CE+0 

.480180E+0 

.275070E+0 

.433840E+0 
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TABLE  VI -D 


PRDBLEM    RECRUIT       «SAMPLE« 


MINIMUM    COST    OF    OBJECTIVE    FUNCTION    IS         . 690549E+04 


GREAT    LAKES    PERCENT 


.330E+02 


BASIS  VECTORS  AND  COEFFICIENTS 


VECTDR 


COEFFICIENT  (X-ZERO  COMPONENT) 


P(  1) 
P{  3) 
P(U7) 
P(  7) 
P(  9) 
P(  12) 
P  (  1  3) 
P(15) 
P(  17) 
P(20) 
PU3) 
P(24) 
P(37) 
P(28) 
(30) 
(32) 
(26) 
(36) 
(22) 
(40) 
PU2) 
P(  6) 
P(U6) 
P(  18) 
P(50) 
P(52) 
(54) 
(5  6) 
(58) 
(60) 
(62) 
P(6U) 
P(66) 
P(68) 
P(70) 
P(72) 
P(7U) 
P(76) 
P(78) 
P(80) 
P(33) 


.256060E+0 

.4460U0E+0 

. 192276E+0 

.735140E+0 

.  136107E+0 

.139941E+0 

. 1399U1E+0 

.  139941E+0 

.39U550E+00 

.205UU3E+0 

.U338U0E+0 

.2UU575E+0 

.488250E+0 

.205UU3E+0 

.147832E+0 

. 167400E+0 

.135875E+0 

. 139500E+0 

•1U7832E+0 

.306900E+0 

.456274E+0 

.21U760E+0 

.  16633UE  +  0 

.436150E+00 

.277732E+0 

.180803E+0 

.186418E+0 

.250429E+0 

.338023E+0 

. 167327E+0 

. 1U3260E+0 

.188708E+0 

.253916E+0 

. 13634UE+0 

.12UU88E+0 

.141284E+0 

.620100E+0 

.214650E  +  0' 

.U8018CE+0 

.275070E+0 

.292950E+0 
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TABIE  VI -E 


PRDBLEM    RECRUIT       "SAMPLE* 


MINIMUM    COST    OF    OBJECTIVE    FUNCTION    IS         .680941E+04 


GREAT    LAKES    PERCENT 


.340E+02 


BASIS  VECTORS  AND  COEFFICIENTS 


VECTDR 


COEFFICIENT  (X-ZERO  COMPONENT) 


(  1) 
(  3) 

(Ul) 
(  7) 
(  9) 
(47) 
(  13) 
(15) 
(17) 
(20) 
(22) 
(2  4) 
(37) 
(28) 
(4  3) 
(32) 
(26) 
(36) 
(30) 
(40) 
(42) 
(  12) 
(4  6) 
(  6) 
(50) 
(52) 
(54) 
(56) 
(58) 
(60) 
(62) 
(64) 
(66) 
(68) 
(70) 
(72) 
(74) 
(76) 
(78) 
(80) 
(33) 


.256060 
.4460-40 
.563850 
.735140 
.  136107 
.192276 
. 139941 
. 139941 
.830700 
.205443 
. 147832 
.244575 
.488250 
.205443 
.433840 
.  167400 
.135875 
.139500 
.  147832 
.306900 
.399889 
. 139941 
.  166334 
.214760 
.277732 
. 180803 
. 186418 
.250429 
.338023 
.167327 
.  143260 
.188708 
.253916 
.136344 
. 124488 
. 141284 
.620100 
.214650 
.480130 
.275070 
.292950 


E  +  0 
-  +  0  . 

+  00 

+  0* 

+  0 

+  0 
0 
0. 
00 
0' 
0 
0 


E 
E 
E 

E+ 
E  + 
E  + 
E  + 
E+0 
E  +  0 
E  +  0 
E  +  0 
E+0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E+0 
E  +  0 
E  +  0 
E  +  0 
E+0 
E  +  0 
E  +  0 
E+0 
E  +  0 
E  +  0 
E  +  0 
E  +  0 
E+0 
E  +  0 
E  +  0 
E+0 
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TABLE  VI -F 


PRDBLEM    RECRUIT       »SAMPLE« 


MINIMUM    COST    OF    OBJECTIVE    FUNCTION    IS         .671361E+04 


GREAT    LAKES    PERCENT 


.350E+02 


BASIS    VECTORS    AND    COEFFICIENTS 


VECT3R 


COEFFICIENT  (X-ZERO  COMPONENT) 


(  1) 

(  3) 

(41) 

(  7) 

(  9) 

(U7) 

(  13) 

(  15) 

(17) 

P(37) 

P(43) 

P(24) 

P(26) 

P(28) 

P(30) 

(32) 

(20) 

(36) 

(22) 

(40) 

(42) 

P(  12) 

P(46) 

P(  6) 

P(50) 

P(52) 

P(54) 

P(56) 

(58) 

(60) 

(62) 

(6U) 

(66) 

(68) 

(70) 

(72) 

(74) 

P(76) 

P(78) 

P(80) 

P(33) 


.256060E+0 

.4U6040E+0 

.  156385E  +  0 

.735140E+0 

.  136107E+0 

. 192276E+0 

.139941E+0 

. 139941 E+0 

.83070GE+00 

.486250E+0 

.433840E+0 

.244575E+0 

.135875E+0 

.205443E+0 

.  147832E+0 

.  167U00E  +  0 

.205443E+0 

. 139500E+0 

.  147832E  +  0 

.306900E+0 

.299889E+0 

. 139941E+0 

.  16633UE  +  0 

.214760E+0 

.277732E+0 

. 180803E+0 

.  186418E+0 

.250429E+0 

.338023E+0 

. 167327E+0 

.145260E+0 

. 188708E+0 

.255916E+0 

.  136344E  +  0 

. 124488E+0 

. 141284E+0 

.620100E+0 

.214650E+0 

.430180E+0 

.275070E+0 

.292950E+0 
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APPENDIX  A 

BASIC  SYMBOLOGY  OF  THE 
TRANSPORTATION  PROBLEM 

It  is  desired  to  ship  a  homogeneous  product  in  integral  amounts, 

aj  (i  =  1,...m)  from  each  of  m  origins  and  to  have  it  arrive  at  n 

destinations  in  such  a  manner  that  each  destination  receives  an  amount, 

b-j  (j  ~  1,...n).  The  cost  Cjj  of  sending  a  unit  amount  from  the  i^n 

origin  to  the  j   destination  is  known  for  all  ij  combinations.  Here 

the  total  to  be  sent  equals  the  total  demanded,  i.e., 

m        n 

y      y- 

i=1  1     j=1  J 

The  problem  is  to  select  the  amounts  x. .  to  be  shipped  from  the 

•  J 

i       origin  to  the  j       destinations  in  such  a  way  that  the  total  cost  is 
minimized.     There  can  be  no  shipments  from  destinations  to  origins, 
which  is  to  say  x^   -  0  for  all  ij.     Symbolically,  we  may  state  the 
problem 

letting: 


m      n 


=  1    L   Oij  XM 
i=1  3=1    iJ    1J 


where         xi-|  2"°  for  &^  *  and  J 


(A-1) 


nimize 

z  subject 

m 

to 

the 

restraints 

h   xij  = 

i=1 

bJ 

j  =  V 

.  «n 

n 
L   x^  = 

ai 

i  =  1,. 

.  ,m 

(A-2) 
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m 


1=1      5=1 


xij  -° 


for  all  i  and  j 


Restating  the  transportation  problem  as  a  standard  linear  program 
in  the  matrix  algebra  form 

minimize      z  =  CTX  (A-1 ' ) 

subject  to  the  restraint 

where        AX  =  B  (A-2«) 


LCij 


a  transposed  vector  of  c. .'s  (all  vectors  here  are  assumed  to  be  verti- 
cal arrays,  hence,  any  horizontal  array  is  a  transpose) „ 


u 


A  = 


1100  . 
0011  . 


00   .  .  .11 

1010  .  .  101 


(A-3) 


the  matrix  of  coefficients  of  the  system  of  equations: 


I 


xij  "  ai 


i  = 


,xa 


I  , 


—  b„ 
i1    1 
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we  delete  the  equation  for  D2,   since  any  one   of  the  equations  may  be 
considered  redundant,  and 


B     = 


J 


the  adjoined  vectors   of  supplies  and  demands. 

The  investigations  of  sensitivity  can  be  stated  in  the  following 
ways: 
I.     Variations  in  Cx 

Assuming  the  already-found  minimum  basic  feasible  solution  is 
XQ,   since  it  remains  a  minimum  basic  feasible  solution 

if  and  only  if 


;j 


*    CJ 


for  all  j 


where 


o  o 


xj  =  p"1pj 


,-1 


and  P       is  the  inverse   of  P,  the   (m  +  n)  x  (m  +  n)   matrix  composed  of 
column  vectors   of  A. 

Then  it  is  true  that  if 


(AC)T  Xj-  AO,  =  .,  -  0j 


the  solution  X  remains  minimum,  i.e.,  in  investigating  the  effect  of 
cost  on  the  minimum  we  first  solve  (A-1)  to  see  if  our  old  solution  is 
still  good.  If  not,  we  must  resolve. 
II.  Changes  in  B 

Changes  in  B,  principally  in  the  partition  of  it  involving  the 
amounts  to  be  supplied,  are  of  the  most  interest  in  this  problem  since 
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there  is  some  probability  that  individual  recruit  stations  will  not 
meet  their  quotas. 

Since   the  altered  B,  B1   may  be  described     B'   =  B  +  A    B  =  P' 
changes  in  the  minimum  basic  feasible   solution  satisfy 

x«  =x0+  ax0 

Xo  =  ^    Po 

and  X^  will  be  a  basic  feasible  solution,   if  and   only  if 

xl.  >»  0  for  all  xl  members   of  X' 

J  jo 

The  changed  minimum  cost   (objective  function)      z'      is 

-T 

^0  ~  u     Ao 

=  zQ  +  CT    A  Xo 
i.e.,  A  *0  =  C    A    X0  (A-3) 

Since  the  values  of  B  are  subject  to  statistical  variation,  the 
planner  responsible  for  budgeting  transportation  will  want  to  know  how 
such  variation  will  affect  his  cost,  and  equation  (A-3)  gives  an  indi- 
cation of  how  to  go  about  it,  as  we  observe  later  in  Chapter  U9 
Statistical  Variation.  For  a  complete  discussion  of  sensitivity  analysis, 
see  Garvin  3  . 
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APPENDIX  B 
PROBABILISTIC  FORMULATION 

In  this  discussion,  extensive  use  will  be  made  of  the  symbols  and 
ideas  presented  in  Appendix  A  and  some  new  ones  will  be  introduced.  A 
full  list  of  symbol  definitions  is  presented  in  the  Glossary,  Appendix 
B. 

In  the  problem  under  consideration  we  are  aided  in  our  probability 
investigations  by  the  fact  that  of  the  system  of  equations  describing 
the  transportation  problem,  there  is  always  one  redundant  equation 
which  may  be  discarded.  Since  the  inputs  from  the  RS's  are  the  members 
of  B  of  most  interest,  we  discard  an  equation  which  describes  the  input 
to  one  of  the  RTC's.  We  now  assume  that  the  distributions  of  the  RS 
outputs,  F~(b),  are  independent  and  of  the  same  form  although  they 
have  different  parameters.  Of  course,  the  assumption  of  independence 
may  be  attacked,  as  it  implies  that  the  events  being  assigned  the 
numerical  values  b,  i.e.,  the  number  of  enlistments  at  the  RS's,  do  not 
depend  on  any  common  phenomenon  such  as  prosperity,  war  fever,  etc. 
(We  do  not  include  the  effect  of  seasonal  fluctuations  since,  although 
they  probably  affect  the  rate  of  enlistment  over  time  intervals  shorter 
than  a  year,  they  probably  have  no  effect  on  the  annual  rate). 

A  complication  of  notation  is  that  when  we  formulate  the  problem 
in  the  simplex  notation,  we  note  that 


B  = 


yS 


corresponds  to 


am 
Lbd 


of  the 


transportation  notation  and  the  b.    of  the  latter  are  restricted  so  that 

J 


y    y 

l—  Q.A    —    L-j  b,« 


J 

.] 
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or  in  the   notation  of   our  probabilistic  simplex 

m  m+ri 


1=1  .i=m+1 


Fortunately  in  this  problem  (as  indicated  above)    n  =  2  and  since   the 

th 
(m  +  n)       equation  is  redundant  and  is  discarded,   there  is   only   one 

random  variable    /i>,-,    namely  [sm+'\ »   tb.at  is  a  function  of  the   others. 

The  vector  B  becomes 


[Al 

rA      " 

& 

A 

A 

= 

A 

0* 

A 

P^ 

m 

rl  A 

3=1 

Since  we   know  the  proportion   /  ,    of  the  total  number   of  recruits    L^jO^ 
which  the  first  Recruit  Training  Center  is  assigned,  we  have  written 


to 


/ti  =  r£A  • 


j=i 

In  order  to  obtain  the  form  of  the  probability  distribution  func- 
tion for  Lj.f   j  =  1,...m  the  following  reasoning  was  used; 

A  recruiter  has  some  probability,  p,  of  enlisting  any  candidate 
he  interviews.  If  these  probabilities  are  assumed  identical  for  any 
one  Recruiting  Station,  the  number  enlisted  in  N  interviews  is  a  ran- 
dom variable  described  by  a  binomial  distribution  with  parameters  N 
and  p.  However,  the  number  of  interviews  held  is  largely  out  of  the 


& 


control  of  the  station,   that  is,   it  is   controlled  in  only  a  general 
way  and  is   not  predictable.     Energetic  RS  personnel  can  raise  the   local 
level  of  interest  but  must  wait  for  candidates  to  show  that  interest. 
If  time  periods  are  properly  chosen  it  may  be  said  that  the  number  of 
people  arriving  to  be  interviewed  at  an  RS  is  a  random  variable   obey- 
ing the  following  restrictions! 

1.  The  probability  of  one  arrival  in  a  period  At  in  length  is 
roughly  proportional  to  the  average  rate  of  arrivals  and  to 
the   length  of  At. 

2.  The  probability  of  zero  arrivals  in  a  period  At  long  is 
roughly  one  minus  the  probability  for  one  arrival. 

3.  The  probability  of  two  or  more  arrivals  in  a  period  At  is 
of  a  much  smaller  order  than  those  in  1   and  2. 

U.  The  arrival  of  any  one  person  is  independent  of  that  of  any 
other. 

These  are  the  standard  assumptions  for  the  establishing  of  a 
Poisson  probability  distribution  law,  and  the  steps  used  in  arriving 
at  it  are  discussed  in  elementary  probability  texts   (vide   Parzen  [5_j). 
We  will  simply  say  that  N   (the  number   of  interviews  held  in  a  RS  in 
some   convenient  time   interval  A  t)   is  described  by  the  mass  function 

-at  (at)n 


P  [N=  n]   *■■"**  MSL  n  =  0,1,2,... 


n   ! 

=  0  otherwise 

where  OC   is  the  average  number  of  arrivals  per  period  (e.g.,  per 
recruiting  day) . 

Now  the  number  of  men  actually  enlisted  in  N  interviews  is  seen 
to  be  a  random  variable  described  by  a  conditional  distribution  function, 
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For  the  j™  station 


(bJm)  =  p 


7  P. In. 


A  -  bj 


N  =  n 


There  is  one  more  complicating  detail,  which  is  that  the  re- 
cruiter is  not  able  to  stop  recruiting  until  he  has  filled  his  quota, 
q..   In  effect,  this  truncates  the  binomial  distribution  governing  the 

J 

"■output"  of  his  n.  interviews,  so  that 


and 


if  n3  "=  qj 

if  q  .  -*  n 


J 


Dropping  the  j   subscripts  for  clarity,   then 

FB|N   (q|n)   =    L  ,    P0>|n) 

BIN  (b|p(b|n)^0} 


ll(f)^1 


i=0 


if  n^q 


h  K, 


Z  ( ?  i  p1^ 

i=0  »     ' 


if  q<  n 


so  that 


Kh    =  1 


and 


where 


K~  = 


1 


£|?]pV-i 


r     = 


i=0 

(1   -  p) 


U 


now  Y0    (b)  =      /_,  F/Q|N   (bln)   V^) 

P  {n]p(n)-0}PIN  N 


=  L  (cxt)nV-  2- (J *r 

n=0  nI       i=0  l     ' 


oo 


-at 


+  I  (at)ni" 


n=q+1 


I(°)pir»-i 


n!  q 

I 

i=0 


Zjjjpim-i 


=   1 
i=0 


oo 


I  (0(t)<:i|pV"i  +  Hii 


I    (O^V^pir^ 


n=0 


q 

k=0 


IUIpV- 


The  size  and  complexity  of  this  expression  suggests  the  desirability 
of  resorting  to  a  computer  when  computations  must  be  carried  out,  and 
for  some  of  our  purposes  we  can  do  so  by  use  of  Monte  Carlo  simulation 
techniques. 

The  planner  responsible  for  budgeting  the  transportation  of  re- 
cruits will  wish  to  know  not  only  the  "best"  way  to  transport  them,  but 
also  what  the  effort  is  going  to  cost.  A  common  measure  of  cost  used 
in  such  situations  as  this,  where  there  is  uncertainty  involved,  is  to 
determine  the  expected  value  of  the  cost;  in  this  case  it  would  be  of 
minimum  cost.  Ideally  the  entire  probability  distribution  of  minimum 
cost  is  desired,  but  problems  of  this  nature  are  often  found  to  be 
too  difficult  to  solve  entirely.  In  this  case  even  if  we  could  solve 
for  the  expected  value  of  minimum  cost,  we  would  not  have  arrived  at 
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a  very  satisfactory  answer  to  the  question,  "How  cheaply  can  this 
transportation  be  carried  out?"  The  reason  the  answer  is  unsatisfac- 
tory is  that  implicit  in  the  use  of  the  expected  value  of  minimum 
cost  to  answer  the  question,  is  the  assumption  that  a  real  system 
could  be  set  up  to  handle  any  transportation  input  optimally.  In- 
stead  of  taking  such  an  ideallistic  approach  we  have  decided  to  be 
satisfied  with  trying  to  find  a  solution  which  will  deviate  least  from 
this  ideal.  The  "analytical1*  expression  for  the  probability  distribu- 
tion of  recruiting  station  inputs  to  the  transportation  problem  (de- 
rived above)  gives  some  idea  of  the  complexity  of  an  analytical  approach 
to  finding  this  expected  value.  In  order  to  avoid  the  difficulty,  the 
FORTRAN  program  for  solving  transportation  problems  was  modified  to 
solve  a  large  number  of  randomised  problems  when  given  inputs  of  quotas, 
probability  of  interview  success,  and  average  rate  of  arrivals  for  each 
RS.  Thus,  it  was  possible  to  generate  a  vector  of  random  inputs  and  to 
obtain  a  sample  mean  value  of  minimum  costs.  In  order  to  give  an  in- 
dication of  how  the  distribution  of  minimum  costs  varies  the  program 
was  designed  to  solve  for  sample  deviation  also.  Here  again  it  should 
be  noted  that  since  many  data  were  picked  at  random  (i.e.,  the  values 
of  p.,  the  probability  of  interview  success,  and  OC.  the  average  arrival 
rate  for  the  j*""  RS),  the  solution  is  only  an  illustration  of  method. 
From  the  previously  mentioned  relation  that 


e[z0(b)]  -s  Zq  (E[bJ)  U.D 


we  see  that  a  solution  of  the  transportation  problem  using  mean  values 
of  inputs  will  give  us  an  upper  bound  on  the  set  of  possible  means  of 
optimal  solutions.  It  is  also  one  specific  description  of  the  input 
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situation  for  which  we  can  find  a  minimum  cost  transportation  pattern, 
but  it  may  not  be  the  "best". 

What  we  seek  for  the  "best"  is  first s  a  loss  function  describing 
the  penalty  we  pay  for  improper  action,  and  second,  a  rule  for  deciding 
what  action  to  take  so  that  we  can  expect  to  pay  the  least  penalty. 
In  other  words,  we  are  trying  to  solve  a  standard  statistical  decision 
problem  and  are  trying  to  determine  an  estimator  Desirable  proper- 
ties of  such  estimators  are  discussed  in  [6j  and  other  standard  sta- 
tistics  texts.  Those  properties  sometimes  conflict  and  are  often 
chosen  somewhat  arbitrarily.   Our  estimator  will  be  chosen  arbitrarily 
also. 

Where  we  first  define 

z  =  min  z(B) 

A 

for  B,  an  estimator  of  B,  a  natural  loss  function,,  „,  would  seem  to  be 

A(B;B)   =  z(BjB)    -  z^    vB1 


where  we  write  z(BjB)  to  emphasize  the  parametric  role  played  by  B. 
Later  it  will  also  be  written  without  the  B. 

The  risk,  or  expected  loss,  we  wish  to  minimize  is 


R(E*jB)  =  e[A(B;B)] 


where  B  represents  the  decision  rule  for  arriving  at  the  final  action, 
the  appropriate  transformation  on  B„ 

A 

Our  rule,  "use  the  estimator  B",  is  that  which  consists  of  acting 
as  though  B  were  the  actual  value  of  B.  The  action  then  is 


a 


z  =  CT  X 


*  =  CT  £~1  B 


where  we  now  regard  z  as  a  transformation  on  B.  The  vector  X,  the 
"minimum  basic  feasible  solution"  is  one  of  the  results  of  the  solution 
of  our  transportation  problem  under  the  restrictions  represented  by  B. 
Equation  (4-.1)  becomes 


min  R  =  min  E  [  z(B)  -  z0  (B)] 


B 

In  words,  equation  (4-.1)  regards  the  loss  as  the  difference  between 
a  fixed  solution  and  an  ideal  solution.  The  fixed  solution  is  that 

A 

arrived  at  by  use   of  the  estimator  B.     The  ideal  solution  is  that  arrived 
at  by  use   of  the  actual  value   of  B. 

Drawing  on  the   notation  of  simplex  algorithm  one  is  tempted  to  say 


min  R  =  min    /cT£_1  (e[b]  )   -  E  [ctP"1b]| 


making  use   of  the   linearity  of  the  transformation  represented  by  P     , 
that  transformation  arrived  at  by  a  simplex  solution  using  the  estima- 
tor  B,     In  fact,   the  authors  took  this  approach  until  certain  computer 
outputs  brought  the  inherent  error   (discussed  below)   to  their  attention. 
Continuing  the  above  reasoning  would  have   led  to  the  following  satisfy- 
ing solution. 

Letting  e[bJ  =  B,  and  substituting,   equation   (a)   becomes 

e  [ctp"1b]  5  C7?"1! 

where   P       is  the  transformation  arrived  at  by  use   of  E[BJin  the  simplex 
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algorithm,  and  by  the  definition  of 


*o< 


B)  =  min  ,fz(B)]  =  min  fcTP"1BJ 
>■         J        p-1 


T— — 1—  T    —  1— 

C   P     B  ^  C   P     B 


-1   j.  -B-1 


for  any  other  P       j^  P     ,   therefore, 


■  min  ?CTP"1I 

A 

B 


[cV1e]] 


min  R     =  [cTP  1B  -  E  [cTP"1b]  j 


A 


i.e.,  B     =     B 


and  P  '  =  P  ' 


the  best  policy  is  to  use  the  mean  value  of  B  as  an  estimator  and  the 
transportation  pattern  computed  from  it. 

The  above  approach  is  in  error  in  that  we  assumed  we  could  use  the 
transformation  r   over  the  entire  space  of  B  vectors  with  some  of  the 
solutions  becoming  non-optimal,  ignoring  the  possibility  that  they 
might  become  infeasible.  However,  the  fact  is  that  the  solution  to 
the  linear  programming  problem  is  either  optimal  or  else  it  is  infeasible, 
where  the  term  feasibility  is  technically  defined.  For  a  solution  to 
be  feasible  it  must  satisfy  the  restrictions 

AX  =  B 
and         x±   5>  0         for  all  i.  U.2) 

For  the  mathematical  details  of  this  aspect  of  the  simplex  algorithm, 
the  reader  is  referred  to  Garvin  [_3J  . 
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We  can  restate  the  definition  of  feasibility  by  defining  the 

—1  —1 

elements  of  the  P~  matrix  as  p.^  and  writing  for  equation  (4.. 2)  (in 

the  case  where  we  started  with  B  and  its  implicit  transformation  P  ) 


m 

I 

k=1 


£A-°  , 


and  we  want  to  know  the  probability  of  this   occurring,   i.e., 

r  m  ~i 

_k=1 

The  implication  of  the  fact  that  P   E  is  optimal  or  else  not 
feasible  is  that  there  are  two  regions  in  the  space  of  B  vectors  and 
they  are  disjoint.  In  one  region  the  loss  function  is  zero-valued 
and  P   is  optimal.  In  the  other,  the  loss  function  takes  on  a  posi- 
tive  value,  so  apparently  the  estimator  B  that  we  seek  is  that  one 
which  is  associated  with  the  region  of  infeasibility  with  the  least 
expectation.  That  is,  the  risk  can  be  more  accurately  defined  as 


R 


(B;B)  =  0  +  J  A(B;B) 
(b  P~1B  (  KJ 


dFr 


where  K  is  the  set  of  all  feasible  solutions  to  the  equation  AX  =  B. 
Thus,  the  estimator  we  seek  is  that  which  minimizes  the  second  term  on 
the  right,  i.e., 

min  R(BjB)  =  min  J    /  A(B;B)dFE  I 
B     P"VKJ       J 

A_-| 

Now,  defining  X*  =  P     B,  and  x*  as  some  vector  value, 
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A  (B;B)  =  0 

>  0 

>[x*£  k]  =  r, 

K 
{X*/K]   =  J  dFx*(x*) 


if  X*£    K 

if  X*£  K,   i.e.,  X*  E  K 


dFx*(x*) 


K 


=  J  f     (x*)dx* 


K 


where  K  is  the  compliment  of  K  in  the  space   of  all  P     B, 

and 

fy*(x*)       =  fx 

1  '20e"m 


fx*(x*)       -  fx   n        j   (xr..xm) 


*&'&""&  (b1,b2*--bm)|j(b1--bm)l 


-=1 


Now 


J(b. , . . ,b  )  = 


3xi  d*i 


3.*-!         3br 


3)xm  ?>x 


m 


rn 


=     P 


I/O 


therefore 


m 


fX*(x*)  :    f 


>i'->*4SPiiXi'*"SPBiXi) 


Here  a  difficulty  arises  due   to  the   non-independence   of   H     .   and 
A,   i  =  1 , . .  .m. 
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However,  assuming  the  foregoing  operations  have  been  carried  out, 
we  can  make  use  of  the  central  limit  theorem  to  say  that  the  random 
variables  S3*t   i  =  1,...m,  can  be  said  to  be  approximately  normally 
distributed  since  they  are  large  sums  of  random  variables,  at  least 
over  long  periods  of  time.  Further,  since  /jm+1  is  also  a  large  sum 
of  these  sums,  it  can  be  approximated  also  as  a  normally  distributed 
random  variable  which  is  independent  of  the  other  m.  Consequently,  we 
are  able  to  write  for  these  approximations 


fB(b)  =  f/31(..^m(bi--V 


m 


£  %{V 


i.e.,  there  is  a  density  function  which  can  serve  our  purpose. 
Now 


E  h(B;B)l  =  JA(B(X);B(X*))dFx*  . 
L     J  -jY#lY*r 


A  ,/\  , 


{x*|x*g  k} 


A 


To  carry  this  out  we  must  redefine  z  to  conform  to  reality,  e.g., 

m 


where 


:(b) 

— 

L  vt 

1=1 

X*' 

1 

= 

0 

if  xf<0 

i 

= 

X*   +    X* 

l 

if  x*     <  0  and  i  is   odd 
1-1 

= 

X*   +    X* 

i         i+1 

if  x*     <  0  and  i  is  even 
i+1 

— 

X* 

otherwise 

(This  definition  of  xf     models  the  action  of  an  RS  in  carrying  out  a 
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fixed  policy).     Finally 

m       m 

P 
dxv;..dxm 


*  in         iu  .  in  iu 

-CO-DO   -OO  /  ,      / 


Here  again  we  have  come  to  a  point  where  for  computational  purposes, 
theory  yields  to  simulation  and  the  theory  merely  provides  a  conceptual 
framework  for  what  we  are  trying  to  do.   Of  the  ideas  written  above,  at 
least  the  loss  function  can  be  useful  in  the  simulation. 

In  investigating  the  results  of  the  above  reasoning  by  Monte  Carlo 
simulation,  a  program  was  designed  to  provide  the  randomized  input  B 
vectors  and  solve  for  their  sample  mean  (one  common  estimator).  This 
mean  vector  was  then  used  as  an  input.  The  program  computed  the  results 
of  transforming  the  randomized  B  vectors  (using  the  transformation  found 
in  the  solution  for  the  mean  value  vector)  altering  it  by  using  the  de- 
finition given  above  for  xf     outside  the  region  of  basic  feasibility. 
It  further  computed  the  loss  function  in  each  case  and  the  sample  mean 
loss  or  sample  risk  for  this  estimator.  In  computing  the  random  B 
vectors,  the  same  concept  of  the  events  giving  rise  to  enlistments  was 
used  as  in  the  analytical  description;  however,  because  of  the  large 
numbers  involved,  and  to  speed  computation  time,  a  normal  approximation 
to  the  Poisson  distribution  was  used.  For  this  approximation  see 
Parzen  5 J,  and  for  the  formula  used  to  compute  random  numbers  according 
to  a  normal  distribution  see  Vaa  I  9  I. 

We  are  able  to  find  sample  risks  given  specific  B  estimators  but 
we  do  not  believe  it  is  possible  to  solve  for  the  best.  The  only  one 
for  which  we  have  solved  so  far  is  the  mean  B  vector,  B.  For  examples 
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of  the  results  of  this  simulation  (using  completely  arbitrary  para- 
meters) see  Tables  VI I -A  through  VII-G,  pages  55  to  61 ,  and  for  a 
discussion  of  the  results  see  Section  6„ 
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APPENDIX  C 
OTHER  POSSIBLE  LINES  OF  INVESTIGATION 
1.   When  a  planner  decides  to  seek  a  solution  to  a  problem  by  using 
the  criterion  of  minimum  expected  loss,  he  is  considering  all  events 
(in  our  case  B  vector  values)  which  have  any  positive  probability  of 
occurring.  Suppose,  however,  he  should  prefer  to  make  a  wager  based 
on  the  assumption  that  in  the  finite  period  of  time  his  system  (the 
U.  S.  Navy)  is  in  operation,  that  events  of  extremely  small  probability 
will  not  occur.  By  basing  his  actions  on  such  an  assumption  he  is  bet- 
ting against  nature  that  he  can  save  an  amount  of  money  in  this  finite 
period  of  time,  realizing  that  he  risks  losing  more  than  the  amount  he 
saves,  but  risks  it  with  a  very  low  probability  that  the  (extra)  loss 
will  occur.  In  effect,  in  our  case  he  wishes  to  find  the  solution 
(B  estimator)  for  the  following  objective  functions 


lin  <  Iz  (BjB)  -  z0(B)J  dF] 


S    JA(B|S)( 

where 

9  represents  the  upper  limit  of  the  small  level  of  probability  of 

the  events  he  wishes  to  bet  will  not  occur, 

z  (BjB)  is  the  objective  function  solution  of  the  transportation 

problem  for  B  evaluated  at  B  and  emphasizes  the  parametric  role 

played  by  B. 

dF,  is  the  differential  of  the  probability  distribution  function 

for  A(Bj§).   A(B?B)@  is  the  9th  percentile  of  A(BjB)  under  the 

probability  law  F^. 

That  the  bet  mentioned  above  might  be  a  desirable  one  to  make 
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could  be  determined  in  the  following  way.     Suppose 9  as  is  likely,   that 
the  expected  values   of  the    /j's,    the  components   of   our  random  vector  B, 
occur  near  their  maximum  values ,   the  q's.     The  reason  this  is  likely 
is  that  the  quotas   of  recruiting  stations  would  probably  be  assigned  in 
just  such  a  manner.     Also  suppose  B  is   our  best  estimator  using  expec- 
tation.    Further j,   suppose  a  sensitivity  analysis  showed  that  small 
changes  of  B  with  values  of  H   increasing  should  cause   the  solution  to 
become  non-optimal,   i.e.,   the   loss  for  small  increments   of  H  increases 
rapidly  and  should  not  satisfy  the   constraint  that  x^  ^  0  for  all  Xj_ 
of  the  solution 5   and  that  equal  small  changes   of  B  with  values   of  R 
decreasing  should  not  cause  the  solution  to  become   non-optimal,   i.e., 
for  small  decrements   of  Rt   the  loss  remains  at  zero.     The  probability 
mass  function  of  one   of  such  /j's,   Vrp°)  would  have  a  form  similar  to 
Figure  A. 


Pgtb) 


tl 


J. 


FIGURE  A 


-J3 


From  the  figure  we  see  that  the  probability  that  R   takes  on 
the  exact  value  (R   +  AjQ)  is  greater  than  the  probability  that  it 
takes  on  the  exact  value  (O    -A/^) 


fe  fi   +  4fi>~Pg  (§  -A£> 


63 


where  A  p  is  some  arbitrary  small  change  in  p  c     This  statement  hinges 

on   our  assumption  of  the  level  of  p   relative  to  q  and  the  fact  that 

q  is  large,   i.e.,   that  H   is   on  the   increasing  side   of  the  mode   of 
p«(b)  and  so  are  both   (^Q  +  Ay0)  and    (y§  -  Ay0)o 

Looking  in  the  same  general  way  at  the  loss  function,  we  see   that 

it  will  probably  have  a  form  similar  to  that  of  Figure  B, 

A(BjB) 


Ix 


■+- 


J3       q  ft 

FIGURE  B 
since  there  will  obviously  be   no  loss  at  H  .     It  might   (and  a  sensitivity 
analysis  would  determine  this)    look  like  Figure  C, 

A(B5B) 


Llx 


T    "T-T-m 


I 


J3        «    J3 

FIGURE  C 
so  that,  in  the  latter  case  the  probability  density  function  of  ,\(BjB) 
might  be  similar  to  Figure  D. 
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PL(A(B;B)) 


lill 


JL 


FIGURE  D 


0 


Even  so  simple  a  discussion  as  that  above  makes  it  clear  that, 
should  these  conditions  hold,  for  all  (or  possibly  even  some  smaller 
number)  of  the  components  of  B,  an  estimator  of  B  with  components  to 
the  right  of  the  prospective  /j's  would  be  a  reasonable  choice,  given 
the  willingness  to  make  the  previously  mentioned  bet  against  nature. 
The  choice  of  such  a  new  estimator,  B,  might  result  in  a  loss  function 
such  as  that  in  Figure  E, 

A(B|B) 


It 


FIGURE  E 


'fi 


and  a  probability  mass  function  similar  to  that  of  Figure  F, 
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PL(  A(BjB)) 


I        t-mi  irtnt 1 


FIGURE  F 


J3 


The  small  standard  deviations  resulting  from  the  parameters  we 
used  in  our  simulations,  indicate  that  this  second  approach  would  not 
be  profitable.  However,  were  the  numbers  of  recruits  sufficiently- 
small  to  "spread  out"*  the  probability  mass  functions  more,  an  inves- 
tigation of  this  additional  approach  might  become  profitable. 

The  effect  of  ignoring  the  tail,  that  is,  those  losses  with 

probability  less  than  0,  is  illustrated  in  Figure  G.  In  effect  the 

A  -»  th 

new  B,  B,  can  be  found  by  renormalizing  F?  truncated  at  the  0  percen- 
tile and  finding  the  new  mean  of  the  region  of  B  vector  which  produce 
the  truncated  function.  An  approximation  to  this  new  mean  (a  new 
sample  mean)  could  easily  be  found  by  use  of  a  computer. 


ft(A(B;b')) 


EZ 


J L 


FIGURE  G 


ys 
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2.       A  possible  improvement  that  might  be  made  in  the  solution  of 
this  problem  of  the   transportation  probability  model  would  be   obtained 
by  investigating  the  problem's  dynamic  aspect.     A  more  sophisticated 
model  than  ours  might  allow  the  planner  to  make  an  optimal  decision  at 
any  point  in  the  time  period   (one 5,   n9    or  an  rt infinite"   number  of  years) 
based  upon  the  conditions  then  prevailing „     It  is  believed  that  a  multi- 
staged  linear  program  such  as  is  discussed  by  Dantzig   U.     and     7  ,    or 
a  dynamic  program  approach  such  as  those  discussed  by  Bellman     8     could 
be   used  to  provide  this  dynamic  aspect. 
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APPENDIX  D 
GLOSSARY  OF  SYMBOLS 


A        The  matrix  of  coefficients  of  the  system  of  linear  equations 
represented  by  equations  (A=2)  and  (A-2').  The  matrix  is 
defined  in  (A-3). 

a        The  total  number  of  recruits  to  be  shipped  from  the  ith  RS, 
of  which  there  are  me 

(X  -  The  average  number  of  candidates  arriving  per  day  at  the 

J       jth  RS„ 

B        The  adjoint  column  vector  composed  of  the  a. '  s  and  b.'s0 

■*■       .1 

B        A  fixed  value  estimator  of  B, 

B        The  vector  of  expected  values  of  the  random  variables  J3a» 

B        A  (second)  fixed  value  estimator  of  B. 

JO-?      The  components  of  the  vector  B  when  the  vector  is  considered 
as  a  vector  of  random  variables 0 

b,       A  fixed  value  in  the  domain  of  JU-j.  Also  the  total  number 
J       of  recruits  to  be  shipped  to  the  yh   RTCS  of  which  there 
are  ne 

CT       The  transposed  column  vector  of  c,  Js, 

C,.      The  total  cost  of  shipping  one  recruit  from  the  i^*1  RS  to 
2  the  j^RTC. 


The  proportion  of  total  input  assigned  to  the  RTC  of  interest, 

E  V      The  expected  value  of  some  random  variable  V  in  which  V  might 
be  a  scalar  or  vector.  In  the  latter  cases  the  symbol  re- 
presents a  vector  of  expected  values „ 

Fv       The  probability  distribution  function  of  some  random 
variable  V, 

J(.eo)    The  Jacobian  determinant  of  a  transformation  with  respect  to 
the  variables  (coo). 

K        The  set  of  all  feasible  solutions  to  the  linear  program  for 
AX  =  B\ 

K        The  complement  of  K  in  the  space  of  all  possible  P   B 
vectors o 
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A(B;B)       A  loss  function^   the  difference  between  the   least  cost 

possible  if  the  problem  could  be  solved   optimally  for  B9 
and  the   cost   obtained  using  the   transformation  derived  by 
using  B  as  the  estimator. 

N.  The   number   of  interviews   occurring  at  the   j       RS. 

-1 

P        The  inverse  of  the  matrix  composed  of  a  set  of  linearly 

independent  vectors ,  used  as  a  transformation  on  B  to 
yield  X. 

P  ,  P    Those  transformations  arrived  at  by  use  of  B  and  B  respec- 
tively in  simplex  solutions  of  the  transportation  problem. 

p..       The  i,k   element  of  matrix  P  . 

P  (v)     The  probability  mass  function  of  some  random  variable  V. 

p.       The  probability  that  the  2       recruiter  will  enlist  a  candi- 
J  date  at  an  interview. 

Q        The  vector  of  q,»'s. 

th 
q.       The  assigned  quota  of  recruits  for  the  j   RS. 
j 

R(BjB)    A  regret  function  or  expected  value  of  loss. 

r  .       r .  =  1  -  p. 

J        J       J 

9  The   level   of  probability  of  the  events  which  the  planner 

bets  will  not  occur. 

X  The   column  vector  composed  of  the  x..'s, 

^  J 

th 

X11  The  total  number  of  recruits  shipped  from  the  i  u  RS  to 

the  jth  RTC. 

XQ       The  vector  of  34-j's  which  yields  the  minimum  z  for  the 
general  case. 

x.  _,      An  element  of  X»„ 
10  o 

X        The  vector  of  Xj_,°'s  which  yields  the  minimum  1   when  B  is 
the  vector  of  restrictions. 

X*       The  vector  valued  random., variable  resulting  from  the  operation 
of  the  transformation  F   on  the  random  vector  B. 

x*       A  possible  value  of  X*„ 
th 

X*       The  i   component  of  X*. 

x*       The  ith  component  of  x*. 
i 
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x?       The  transformed  xf"  which  models  the  actual  situation  of  an 
i  ■*- 

RS  carrying  out  a  policy. 

z        The  objective  function  expressed  in  equations  (A-1 )  and 
(A-1«). 

z        The  transformation  which  yields  the  minimum  objective 
function^  i.e0,  C  Xc 

z0(B)     The  minimizing  transformation  found  using  B  as  a  set  of 

parameters  and  considering  B  as  a  vector  of  random  variables 

z  (B)     The  minimizing  transformation  found  by  using  K  as  the  speci- 
fic value  of  Bo 

z  (B)     The  transformation  of  B  which  gives  the  minimum  of  the 
objective  function  when  acting  on  B„ 
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